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CN I Abstract 

. In this paper we study the simplest deformation on a sequence of orthogonal poly- 

Q> I nomials, namely, replacing the original (or reference) weight wq{x) (supported on M 

' or subsets of M) by wo{x)e~*^. It is a well-known fact that under such a deforma- 

■ tion the recurrence coefficients denoted as and evolve in t according to the 

^ , Toda equations, giving rise to the time dependent orthogonal polynomials, using Sogo's 

?H ' terminology. 
■ - - ' If ujQ is the Gaussian density e~^ , x € M, or the Gamma density x°e~^ , x € M+ , 

a > —1 , then the initial value problem of the Toda equations can be trivially solved. 
This is because under elementary scaling and translation the orthogonality relations 
reduce to the original ones. However, if wq is the Beta density (1 — x)"(l + x)^ , 
X G [—1, 1] , a,l3 > —1, the resulting "time-dependent" Jacobi polynomials will again 
satisfy a linear second order ode, but no longer in the Sturm-Louville form. This 
deformation induces an irregular point at infinity in addition to three regular singular 
points of the hyper geometric equation satisfied by the Jacobi polynomials. 

We will show that the coefficients of this ode are intimately related to a particular 
Painleve V. In particular we show that Pi{n,t) , where Pi{n,t) is the coefficient of 
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of the monic orthogonal polynomials associated with the "time-dependent" Jacobi 
weight, satisfies, up to a translation in t, the Jimbo-Miwa cj-form of the same Py; 
while a recurrence coefficient an{t), is up to a translation in t and a linear fractional 
transformation Py (q:^/2, — /3^/2, 2n + 1 + a + /3, —1/2). These results are found from 
combining a pair of non-linear difference equations and a pair of Toda equations. 

This will in turn allow us to show that a certain Fredholm determinant related to 
a class of Toeplitz plus Hankel operators has a connection to a Painleve equation. 

The case with a = /? = — 1/2 arose from a certain integrable system and this was 
brought to our attention by A. P. Veselov. 



j,k=o 



1 Introduction 

The study of Hankel determinants has seen a flurry of activity in recent years in part due to 
connections with random matrix theory. In this paper we will investigate Hankel determi- 
nants for a weight of the form 

(l-x)"(l + x)^e-*^ 

on the interval [—1,1]. Here we take i e IR. We will call this a time dependent Jacobi 
weight. Our ultimate goal is to produce a non-linear second order differential equation that 
is satisfied by the logarithmic derivative of Dn{t) , where Dn{t) is the determinant of the 
Hankel matrix generated from the moments of the weight: 

Dr,{t) := det := dct (^j' x^+^ (1 - x)°(l + xfe-'^d:^ 

and we shall initially assume a, /9 > 0. 

The moments ^Xkif) can be evaluated as follows: 

lik{t) = {-lf —^lioit), /c = 0,l,2,... 

where 

^o(t) = 2"+^+^r(a + l)V{p + 1) e* M(/5 + 1; a + /3 + 2; -2t), 
and M(a; 6; z) is the Kummcr function with parameters a and h. Because 

(a) 

^M(a; h- z) = j^^M{a + k;b + k;z), 

we find, 

l,,{t) = 2"+^+^r(a + l)r(/3 + 1) e^ ^ (^^^ (-2)^ ^ + m(/3 + r + 1; a + /j + r + 2; 2t) . 
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Our interest in evaluating Dn[t) comes from the fact Dn{t)/Dn{0) is the generating function 
of the hnear statistics 

n 

and can also be thought of as the partition function for the random matrix ensemble with 
eigenvalue distribution 

n 

n - ^'^)' 11(1 - xini + xif e-'^^dxi. 

l<j<k<n 1=1 

The paper is divided into the sections as follows. In the next section we reproduce known 
results that we call coupled Toda equations. We include them in this paper as reference. In 
section 3 we consider ladder operators and derive fundamental equations that are the basis 
for everything that follows. They give us coupled difference equations and coupled Riccati 
equations in certain auxiliary quantities denoted as r„(t) and Rn{t)- 

The Ricatti equations allow us to find a non-linear second order differential equation that 
is satisfied by the recurrence coefficient Q;„(i). A rational change of variable applied to CKn{t) 
is then a solution to a Painleve V in standard form and can be found in section 4. 

In section 5 we identity the function which satisfies the continuous and discrete cr -form 
of our Py. We show that the cr -function of Jimbo, Miwa and Okamoto is given by 

t n 
(j{t)^-p^{n,t/2)--t + n{n + P), 

and since ^ 

Pi{n,t) = — log D^{t), 

consequently Dn{t) is related in a simple way to the r -function of the Py. 

In section 6 we show how the Hankcl determinants can be also expressed as determinants 
of finite Toeplitz plus Hankel matrices. For some special cases, these latter determinants 
are known exactly and hence so arc our Hankel determinants. Thus we, in a round about 
way, produce second order equations that have solutions that are logarithmic derivatives 
of Predholm determinants. This should not come as a great surprise as this is a common 
occurrence in random matrix theory for the classical ensembles. 

Finally using known results for the Predholm determinants we are able, in some special 
cases, to write down asymptotic expansions for these determinants and make some predic- 
tions about higher order terms. 

2 Preliminaries: Notations and time evolution 

The purpose of this section is to derive two coupled Toda equations that involve the recur- 
sion coefficients of the time-dependent Jacobi polynomials. This is not a new result. The 
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rather more general To da-hierarchy, can be found for example, in [30], [23] and [37] . Ours 
corresponds to the first of the hierarchy. See [23] for a discussion of this in relation to Sato's 
theory. See also [T| for the "multi-time" approach to matrix models. 

We include the necessary computations here for completeness sake and to set the nota- 
tions to be used throughout this paper. 

To begin we consider general orthogonal polynomials Pi{x) with respect to the weight 
Wo{x)e~^^ on [—1,1]. The weight Wq will be known as the "reference" weight. The orthog- 
onality condition is 



1 

tx 



Piix)Pj{x)woix)e-''^dx = hi{t)6ij, (2.1) 



-1 

and the t dependence through e~*^, induces t dependence on the coefficients. We normalize 
our monic polynomials as 

P^{z) =z^ + pi(n, t)z"-i + ... + P„(0), (2.2) 

although sometime we do not display the t dependence of coefficients of z^~^ . 

An immediate consequence of the orthogonality condition is the three terms recurrence 
relation 

zPn{z) = Pn+l{z) + anPn{z) + /3„P„-l(z) (2.3) 

with the initial conditions 

Po(^) = l, /3oP-i(^) = 0. (2.4) 
An easy consequence of the recurrence relation is 

a„(t) = pi(n,t)-pi(n + l,t), (2.5) 
and a telescopic sum of the above equation (bearing in mind that Pi(0,t) = 0) leaves 

-5^a,(t) = Pi(n,t). (2.6) 

i=o 

First let us discuss the derivatives of a„ and /?„ with respect to t, as this yields the simplest 
equations, where we keep Wq quite general, as long as the moments 

fii{t) := j x'wo{x)e-*''dx, i = 0,l,... (2.7) 

of all orders exist. Taking a derivative of with respect to t 

= - j wo{x)e~^'' xPl{x)dx = -anhn, (2.8) 



i.e., 

(log/i„)' = (2.9) 
and since = hn/hn-i, we have the first Toda equation, 

= K_i - (2.10) 
We define Dn{t) to be the Hankel determinant 

D„(t) = det(/i.+,(t))i;-io. (2.11) 



It is well-known that Dn{t) = YYi=o hi{t). This yields in view of (2.9) that 



Also, 



d """^ 
-logZ}„(t) = -5^a,(t) = Pi(n,t). 

j=0 



d 

= ^ y PnPn^iWoe ^""dx 



xPnPn-iWoe ^""dx + K-i — Pi(ra,t) 
-K + — Pi(ra,t), 



and therefore 



^Pi(n,t)=/3„(t). (2.12) 



But since cxn — Pi(^) — Pi('^ "I" 1); we have the second Toda equation, 

a'^ = (3n- (3n+i. (2.13) 
To summarize we have the following theorem. 

Theorem 1 The recursion coefficients a„(t) and (3n{t) satisfy the coupled Toda equations, 

(3'^ = {an-i-an)l3n, (2.14) 

a'^ = Pn-Pn+i. (2.15) 

It is also worth pointing out that in view of (2.11) we have the obvious Toda molecule 
equation 

n ^ t +\ a u\ Dn+i{t)Dn-i{t) 
^log/^„(t) = j^pM.t)=m = . 
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3 Ladder operators, compatibility conditions, and dif- 
ference equations. 

In this section we give an account for a recursive algorithm for the determination of the 
recurrence coefficients /?„ based a pair of ladder operators and the associated supple- 
mentary conditions. Such operators have been derived by various authors over many years. 
Here we provide a brief guide to the relevant literature, [7], [S], [H], [ID], [H], [13], [15], [20], 
|21j and [26j. In fact Magnus in [26] traced this back to Laguerre. We find the form of the 
ladder operators set out below convenient to use. 

For a sufficiently well-behaved weight (see [15] for a precise statement) of the form w{x) = 
Q-v(x) lowering and raising operators are 

P^(Z) = ~Bn{z)Pn{z)+PnAniz)P^^iiz), (3.1) 

P^_i(^) = [5„(z) + v'(^)]P„_i(^)-A„_i(^)P„(z), (3.2) 



Aniz) := ^ f "^^^^^^Pl{y)w{y)dy, (3.3) 



where 



K J-i z-y 

Bn{z) := ^ "'^""^ Pr.{y)Pr,^,{y)w{y)dy. (3.4) 

K-i J-i z-y 

(3.5) 

Here we have assumed that w{±l) = 0. Additional terms would have to be included in the 
definitions of A„(z) and Bn{z) if w{±l) ^ 0. See [13] and [T5] . 

A direct calculation produces two fundamental supplementary (compatibility) conditions 
valid for all z] 

Bn+iiz) + B^{z) = {z- an)A^{z) - y'{z) (Si) 

1 + {Z- an)iBn+liz) - Bn{z)) = Pri+lAn+l{z) - PrAn-liz). (82) 

We note here that (5*1) and {S2) have been applied to random matrix theory in [36]. It 
turns out that there is an equation which gives better insight into the and if (5*1) 
and (S'2) are suitably combined. See [T7] . 

Multiplying (S'2) by An{z) we see that the r.h.s. of the resulting equation is a first order 
difference, while the l.h.s., with {z — a„) replaced by Bn+i{z) + Bn{z) + v'(z), is a first 
order difference plus An{z) . Taking a telescope sum together with the appropriate "initial 
condition" , 

Bo{z) = A^i{z) = 0, 

produces, 

n-l 

Bl{z) + v'{z)Bn{z) +J2Aj{z) = PnAn{z)An^,{z). {S'^). 

j=0 
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This last equation will be highly useful in what follows. The equations (S'2) and (5*2) 

were also stated to See also [22]. 

If Wq is modified by the multiplication of "singular" factors such as |x — t|" or a + 
bH{x — t), where H is the unit step function, then the ladder operator relations, (S'2) 
and (S'2) remain valid with appropriate adjustments. See [19], [18], and [1]. 

Let ^E'(-z) = Pn{z) ■ Eliminating Pn-i{z) from the raising and lowering operators gives 

where we have used (S'2) to simplify the coefficient of \E' in 03.61) . 

For the problem at hand, w{x) is the "time-dependent" Jacobi weight, and now we must 
suppose that a > and /5 > so that our weight is suitably well-behaved. Then 

w{x) := (1 -x)"(l + x)V*^ xG[-1,1], (3.7) 
^/(z) := -a\og{l - z) - plog{l + z) +tz, 

= 7 + 

z — 1 z + 1 

z~y + + + 

Substituting these into the definitions of An{z) and Bn{z) and integrating by parts, pro- 
duces, 

^""^'^ - "731+ ^ + 1 ' 

BJz) = ^ + ^ 



z + 1 ' 



where 



Rn{t) := ^ C ^ (1 - yni + yfe-'ydy, 
K J-i l-y 



K^i J-i l-y 

Substituting the expressions for An{z) and Bn{z) into (Si) and (S'2) , which are identities 
in z , and equating the residues of the poles at 2; = ±1, we find four distinct difference 
equations and one which importantly performs the summation Ylj Rj • 

- (r„+i + r„) = a- i?„(l - a„) (3.9) 
r„+i + r„ = 2n + l + P-{Rn + t){l + an) (3.10) 
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rl + arn = PnRnRn-l (3.11) 
(r„ - n)2 - /3(r, - n) = /3„(i?„ + t)(i?„_i+t) (3.12) 

f j3 — a\ an 

1^ — J'^""^^ ~ t^n-rn[rn-n)- 2_^Rj (3.13) 

= —^[Rn{Rn~l +t) + {Rn + t)Rn-l\- 

We now manipulate the equations fl3.9l) - fl3.13l) with the aim of expressing the recurrence 
coefficients a„, /5„ in terms of r^, Rn , and of course n, t. 
Adding ([33D and (l3J0D yields, 

2Rn = 2n + a + 13 + l-t-tan, (3.14) 

thus, is "easily" expressed in terms of 
Subtracting (13. 9p from (13.101) gives, 

(3-a + l-t ft ^\ 
r„+i + r„ = nH I - + 1 a„. (3.15) 

Eliminating PnRnRn-i from (13.111) and (I3.12p we find, 

n{n + (3)-{2n + a + (3)rn = /?„[t' + t(/2n-i + Rn)]- (3.16) 
Now with the aid of (3.11), replacing j3nRn-i by (r^ + arn)/Rn, in (3.16) we find, 

t(t + Rn)Pn = n{n + P)-{2n + a + /3)r„ - ^(r^ + ar„), (3.17) 

-fin 

and this expresses /3„ in terms of i?„, r„, n, t . It is important to note the absence of Rn±i, 
and rn±i- The reader will notice that the above manipulations prove that we have expressed 
an and the recurrence coefficients in terms of auxiliary quantities r„ and Rn ■ 
This is summarized in: 

Theorem 2 With , Rn , an and Pn as defined above and with a,P >0 

tan = 2n + l + a + f3 -t-2Rn, (3.18) 

t{t + Rn)Pn = n{n + (3)-{2n + a + (3)rn - ^{rl + ar„), (3.19) 
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In what follows we will find two Riccati equations, one in r„ with coefficients involving Rn 
and another with the roles of r„ and Rn reversed. We will first show that 

d "—1 

-logDjt) =n-2r„-/3„t = = p^{n,t). (3.20) 

j=Q 

The equation (3.20) can be derived as follows: Replace n by j in (13.181) and sum over j 
from to 77,-1 to obtain 

^ ^ nin — 1) nia + (3 + 1 — t) t 

Rn = — + ^ > a,-. 

J 2 2 2 ^ ^ 

j=0 j=0 

Now we can obtain from (I3.13P another expression for Rj, 

f (3 — a\ an ^ tjSn . „ t-, \ 2 

Rj = \— — jrn + — -trn-r^ + nrn + —{Rn + Rn-i)+rn + <yrn 



+ n - t r„ + — + —{Rn + Rn-lj 



j=0 

/ a 

V 2 7 " ■ 2 ■ 2 

= "^"+; + ^^-tr.-^, (3.21) 

where we have ehminated tPn{Rn + Rn-i) using (3.16). These last two equations yield 

n-l 

Pi(n, t) = - ^ a^- = n - 2r„ - 

i=o 

/^From this we also deduce (see (2.5)), 

an = 2(r„+i - r„) + t{(3n+i - /?„) - 1- (3.22) 

Hence, in view of (2.12), 

an = 2(r„+i - rn) - - 1, 



(3.23) 



I.e., 

dan , P — a + 1 — t ft \ 
^~(U + + ^ = 2 I riH I - + i?„ 1 a„ - 2r„ , , 

where we have replaced r„+i by r„ plus an additional term with (3.15). Finally, 

t— |-a„ = 2n + p — a — t — (t + 2it„)Q;„ — 4r„, 

dt 



or replacing q;„ in favor of i?„ from (3.18), 
dR 



t- 



dt 



at + {2n+l + a + P - 2t)Rn - 2Rf^ + 2tr„ 



(3.24) 



This is a Riccati equation in i?„ , but with r„ appearing hnearly. 



ra-l 



Now since 

sec (2.6) and (2.11), we find, upon taking a derivative of (3.20) with respect to t, 

drn „ dPn 



or 



(3.25) 



where we have replaced /5'j(t) by (a„_i — an)/3n with the aid of the first Toda equation, 
(2.13). A simple computation with (3.14) gives, 

t 



:(ttn-l — Oln) — —1 + Rn — Rn-1, 



hence 



— /3n{Rn — Rn-l) 



Pn-R-n 
Rr, 



Rr, 



(3.26) 



t{t + Rn) 



t 



n{n + 13) -{2n + a + l3)rn - — (r^ + ar„) 

Rn 



Rn 



which is a Riccati equation in r„. We summarize in the following theorem: 
Theorem 3 The quantities r„ and Rn satisfy the coupled Riccati equations: 



t- 



dRn 



dt 
drn 



= + (2n + 1 + a + - 2t)Rn - 2Rl + 2r„i, 



dt 



Rn 
t{t + Rn) 



n[n 



+ p)-{2n + a + P)rn - ^ (r^ + ar„) 



rj + Oifn 

Rn 



(3.27) 
(3.28) 



Wc end this section by pointing out that the above equations not only produce differential 
equations in our various unknown quantities, but also a pair of coupled non-hnear first order 
difference equations in Rn and r„. If we substitute (5n into (3.11) we obtain the following 
result. 
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Theorem 4 The quantities r„ and i?„ satisfy the coupled difference equations 

2t(r„+i+r„) = ARl + 2Rr,{2t - 2n - 1 - a - (3) - 2at (3.29) 

n(n + /3)-(2n + a + /5)r„ = (r^ + ar„) f — 4 + ^ + (3-30) 



RnRn — l Rn Rn—1 

together with the "initial" conditions 

ro = 0, (3.31) 

_ a + P + lMil + P;a + P + l;-2t) 
^ ~ 2 M{l + (3;a + (3 + 2;-2ty ^^"^^^ 

The initial condition for Rq can be found by direct integration. Observe that the represen- 
tation of Rq in terms of a ratio of the Kummer functions aUows for the analytic continuation 
of a, P down to a = P = —1/2, due to the relation, 

\imbM{a;b;z) = azM{a + 1,2, z). 

Hence 

lim Ro{t) = -( ^ - 1 

„^_l/2,/3-.-l/2 ' 2 \Io{2t) 

Indeed, by formally continuing (3 so that /5 + 1 = —k, k = 0,1, 2..., we find 

expressed as ratio of Laguerre polynomials of degree k. It is clear that iterating (3.30) and 
(3.31) with the above Rq and ro = will generate rational solutions (in the variable t) 
of our Py derived in section 4. It is interesting to note that Rq is also a rational function 
of a and t, therefore for the values of the parameters stated above our Py is a rational 
function in a and t. 

Also note that the above equations define the quantities r„ and Rn for all a > —1 and 
(3 > —1 . To verify our answers we return to the pure Jacobi case and let t = 0, then (3.14) 
gives, 

, a + (3 + l 
i?„ = nH , 

and is consistent with (3.21) at t = 0. Now equating (3.11) and (3.12) at t = 0, gives, 

n{n + P) 
" a + P + 2n 

and 

+ ar„ 4n(n + a){n + P){n + a + P) 



Pn 



RnRn-1 i2n + a + P + l){2n + a + P-l){2n + a + pf 
With the Rn given above, we find a„ from (3.14) at t = 0, 



{2n + a + p){2n + a + p + 2y 
which are in agreement with those of [T5] . 
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4 Identification of Py 



The idea is to eliminate r„(i) from our coupled Ricatti equations to produce a second order 
ode in Rn{t). This is straightforward and messy and we omit the details. A further change 
of variable, 

Rn{t) - 



t 



leaves, after some simplification. 



where the last term is a particular rational function of two variables defined as follows: 



Rat(|/, t) 
Therefore, 

satisfies. 



2y{y + l) , 2{2n + l + a + (3)y (y-l 



y-1 



+ 



+ 



a 



y{t) -1 + 



Rnity 



^y{y - 1) t t y-i V 



-Try 

2 y 



which is almost a Py. To fit the above to into a Py, we make the replacement t ^ t/2 
followed by 

y{tm = Y{t), 

and find 

t 



2r(r- 1) 

+ {2n + l + a + P) 



2 Y 



t 2 Y -1 



(4.33) 



which is 



Fv{a^/2, -(3^2, 2n + 1 + a + (3,d = -1/2). 
The initial conditions are 

Y(0) = 1, Y'(0) = ^- . 

It is well known that there is a Hamiltonian associated with Py. To identify it, we substitute 



R„{t) := -tq 
rn{t) := -pq{q - 1) + p q, 



(4.34) 
(4.35) 
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where p = pit), q = q(t) into (3.20), and choose p so that the resulting expression is a 
polynomial in p and q. There are two possible p : p = n and p = n + (3. 

Case 1. p = n 

tpi{n, t) + n{n + a + (3) — nt = p{p + 2t)q{q — 1) — 2ntq + (3 pq + a p{q — 1) (4.36) 

Case II. p = n + (3 

tpi{n, t) + n{n + a + 13) — aP — nt 

= p{p + 2t)q{q-l) + 2{p-n)qt + ppq + ap{q-l). (4.37) 

Replacing t by t/2 we see that the l.h.s. of (4.36) and (4.37) are the two Hamiltonian tHi 
and for our Py. The Hamiltonian as presented in Okamoto [31] (see also [32]) is 

tH = p{p + t)q{q - 1) + a2qt - aspq - aip{q - 1), 

where 

a = — , b= — c = ao-a2, " = "2' "0 = 1 - "i - "2 - as- 
Comparing with our tHi, while keeping in mind that the t is in fact t/2, we find 
0^2 = — n, 0^3 = — /3, ai = —a, a;o = ?T, + l + a + /?, 
a = — , h = -'--, c = 2n + l + a + l3. (4.38) 

Comparing with our tH2, we find 

a2 = (3 — n, 0^3 = — /5, ai = —a, ao = n + 1 + a 
a = y, ^ = -Y, c = 2n + l + « + /?. (4.39) 

Hence both Hi and H2 generates our Py, where 

Y{t) = 1 ^ 



g(t/2)' 



5 The continuous and discrete cr-form of Py • 

Recall from the section 2, that 

^logL'„(t) = Pi(n,t), 
13 



and 



p[{n,t) = (34t). 



(5.1) 



Now we come the continuous a -form of Py satisfied by p^(n, t), with n fixed and t being 
the variable. The idea is to express r„, in terms of Pi{n,t) and its derivative with 
respect to t. Let us begin with (3.20), 



Pi(n,t) 



/^From the last equality of (5.2) we have 

rn{t) = - 



n - 2r„ - tl3n 
n-2rn- tp[{n,t). 



n- -^{tPi{n,t)) 



Under some minor rearrangements, equations (13.191) and (13.261) become 

t(3nRn + ^irl + arn) = n{n + (3) - {2n + a + (3)rn - f (3n 



-tf3nRn + ^{rl + ar„) 



tr' 



(5.2) 
(5.3) 

(5.4) 
(5.5) 



Rn we find 



respectively, which is a system of linear equations in l/Rn and Rn- Solving for l/Rn and 
2t 

(r^ + ar„) = tr'^ + n[n + fj) - [2n + a + P)rn - t^'Pn 

(5.6) 



+ ar„) = tr; + n{n + (3) - {2n + a + /3)r„ - t'^pn 

2tpnRn = -K + n{n + f3)-{2n + a + (3)rn-f(3n- 
Taking the product of the above we arrive at an identity free of Rn : 

4t^Pn{rl + arn) = [n{n + p) - {2n + a + /?)r„ - t^Pn? - {Kf. (5.7) 
To identify the a -function of Jimbo and Miwa [25j, we replace t by t/2 so that 

Rnitm = -2(l_V(t))' 
and substitute the above in (3.27) in the variable t/2. After a little simplification we find 

t— = tY- 2r„(t/2)(l - - (r - l){aY + 2n + p). 



Comparing this with the first equation of (C.40) of [25], we have 

-rn(t/2), 
a, 



z{t) 

do — Oi + Goo 



36q + 6i + 6c 



-2n - /3, 
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and consequently 

l-eQ-9i = 2n + l + a + (3 = c, 

consistent with the parameter c of our Vy. Furthermore comparing (4.33) with (C.41) of 
we find a possible identification; 



Oq — 9i + 9c 

a = 

2 

_ 9o — 9i — 9c 
^ ~ 2 



and consequently, 
But since 



9q = — n, 9i = —n — a — P, 9oo = a — j3. 



jMt) = 4t) = -r„(t/2), 



and bearing in mind (5.3) we have, upon integration and fixing a constant, 

1 nt 
a{t) = -tp,{n, t/2) -- + n{n + (3). (5.8) 

The (T-form of our Py is essentially a second order non-linear ode satisfied by p^{n,t), and 
reads 

{tay = [a- ta' + {2n + a + p)a'f + 4[a - n{n + 13)- ta'][{a'f - aa% (5.9) 
with the initial conditions 

a(0)=n(n + /3), ^^(0) = -r„(0) = - ^_[^ . 

a + p + 2n 

After some calculations we find that (5.9) in fact the Jimbo-Miwa a -form (C.45) with 

z/q = 0, z/i = -a, z/2 = n, z/3 = n + /3. (5.10) 
To obtain (5.9) we first replace t by t/2 in (5.7), and substitute 
r„(t/2) = -a'it) 
|r„(t/2) = -a"{t) 

- 2Ap^(.,v2)^4| [-(^)+-^/^^-(-+^)] 

into (5.7) at t/2. Furthermore, since p^(n, t) = (logD„(t))', we have 

^ / N ^ / N r /"* cr(2s) - nin + 3) + ns , " 

D„(t) =D„(0)exp / ^ ^ ds , 

Jo s 
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where -D„(0) given by (1.6) of [3]. 

We expect there exists a discrete analog of the continuous a -form, namely, a difference 
equation in the variable n, satisfied by p^(?T,,t) with t fixed. To simplify notations, we do 
not display the t dependence. The idea is similar to the continuous case namely, we express 
Pn, fn and Rn in terms of Pi(n) and Pi(n ± 1), and substitute these into (3.11) that is, 

+ ar„ = (3nRn Rn-l- 

To begin with, we note that (3.20) is linear in /?„ and r„, which we re-write as; 

t(3n + 2rn = n-p^{n). (5.11) 

We now find a another linear equation in and r„. First note that a„ = Pi{n) — pi{n + l), 
and from (3.14) we have, 

2Rn + t = 2n + l + a + P + t[pi(n + 1) - Pi(n)]. (5.12) 

The sum (5.12) at "n" and the same but at "n — 1", leaves, 

Rn + Rn-l +t=i2n + a + P) + {t/2)[p,{n + 1) - p,{n - 1)]. 

Substituting the above into (3.16) results the other linear equation mentioned above; 

tPn{2n + a + p+ (t/2)[pi(n + 1) - p^{n - 1)]} + {2n + a + /5)r„ = n{n + p). (5.13) 

Solving for tj3n and 2r„ from the linear system (5.12) and (5.13), we find, 

it/2)[n - p^in)][p,in + 1) - p^(n - 1)] - n(n + /3) 



2r. 



tpn 



n+{a + P)/2 + (t/2)[pi(n + 1) - p,(n - 1)] 
[n - p^{n)][n + {a + 13) /2] + n{n + a + (3) 
n+{a + P)/2 + (t/2)[pi(n + 1) - p,(n - 1)] ' 



(5.14) 
(5.15) 



and the discrete cr-form results from substituting (5.12), (5.14) and (5.15) into (3.11). 

Imagine for a moment that we leave our original problem behind and only consider 
the functions Y and a that satisfy the two Painleve equations (4.33) and (5.9) with the 
appropriate initial conditions. Then our orthogonal polynomials Pn{z,t) satisfies the linear 
second order ode 



^"{z) + P{z)^'{z) + Q(^)^(^) = 0, 



where 



P(^) 

Q(^) 



l + a 1 + /3 

\ -t 

z-1 z-1 



z-[i + r(2t)]/[i-r(2t)] 



(5.16) 
(5.17) 



+ 



2{z 
1 



2{z + l 
A 

dt 



n+ia{2t) 



(r{2t) +n{a + l) - 
a{2t) +n{a + l + 2t) 

[1 - r(2t)] 



Y{2t) d 
~2 Jt 



a{2t) 



Id,. 

1 



(1 - Y{2t)) 



z-{l + Y{2t))/{l-Y{2t)) 



(5.18) 
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This is a deformation of the classical ode satisfied by the Jacobi polynomials. When t = 
this ode reduces to a hypergeometric equation. 



6 Toeplitz and Hankel determinants. 

In this section we introduce certain matrices that are combinations of finite Toeplitz and 
Hankel matrices. There are identities that link these matrices directly to the Hankel moment 
matrices that appear in the first section of this paper and define our quantity -D„(t). We 
will use these identities in some special cases to get exact formulas for Dn{t) and, as a 
by-product, find Painleve type results for some other interesting determinants. We include 
the Toeplitz/Hankel computations because as far as we know they are not written down 
explicitly in this form in any other place. However, the Case 2 was established already by 
two of the authors in [Sj Sec. 2]. The current derivation follows that in [5j. 

Given a sequence {ak}'^^_ao '^^ complex numbers, we associate the formal Fourier series 

oo 

a{e'') = "fc^""'' ^ ^- (6-1) 

k=—oo 

The n X n Toeplitz and Hankel matrices with the (Fourier) symbol a are defined by 

Tn{a) = (aj_fc)J;^o, Hn{a) = {aj+k+iY-'klo ■ (6.2) 

Usually a represents an -function defined on the unit circle T = {z G C : |2;| = 1} , in 
which case the numbers are the Fourier coefficients, 

ttk = — [ a{e'^)e-'^Ue, k e Z. (6.3) 

Notice that while the matrices Hn{a) are classically referred to as Hankel matrices they 
are not the same as the Hankel moment matrices considered in the previous sections of this 
paper. To make the connection to Hankel matrices defined by moments we write 

Hn[h] = (&,+fc)"io> hk = - f h{x){2xy+Hx, (6.4) 

where h{x) be an -function defined on [—1,1]. Notice the difference in notation in 
comparison to (12.71) and (12. lip . Our goal in this section is to prove four identities. Let 
z = e^^ . Then for each n ^ 1 the following statements are true. 

1. If a(e^^) =6(cos^)(2 + 2cos0)-i/2(2_2cos^)-i/^ then 

det {Tn{a) - H^iz-^a)^ = det Hn[b]. 
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2. If a{e'^) = b{cos9){2 + 2cos9)-^/'^{2-2cos9y/'^, then 

det (r^ia) + H^{a)^ = det H^[b]. 

3. If a(e^^) =6(cos^)(2 + 2cos^)^/2(2_2cos0)-^/2 ^^en 

det (Tn{a) - iJ„(a)) = det Hn[b]. 

4. If a(e^^) =6(cos0)(2 + 2cose)i/2(2-2cos^)i/2 then 

^ det (T„(a) + H^{za)^ = det 



In these identities the function a is always even, which means in terms of its Fourier co- 
effcients that Ofc = ci-fc . Moreover, in these formulas we assume that a G (T) , which 
implies that (and in Case 4 is equivalent to) b G L^[—l, 1] . We also remark that Cases 2 
and 3 can be derived from each other by making the substitutions a(e*^) i— >■ a(e**-^"'"'^^) and 
b{x) ^— b{—x) . 

Our interest in the above formulas stems from the circumstance that they allow us to 
use existing results [6] on the asymptotics of the Toeplitz+Hankel determinants with well- 
behaved symbols a in order to derive the asymptotics of the Hankel moment determinants. 

In the above identities four types of finite symmetric Toeplitz+Hankel matrices as well as 
a finite Hankel moment matrix occur. These finite matrices can be obtained from their (on- 
sided) infinite matrix versions by taking the finite sections. It turns out that these infinite 
matrices are related to each other in a very simply way, namely they can be transformed into 
one another by multiplying with appropriate upper and lower triangular (infinite) matrices 
form the left and right. These identities for the infinite matrices will be established in the 
next theorem (and the remarks afterwards) in most general setting, where we do not assume 
that the symbols are -functions. 

Let us introduce the infinite matrices 



/ 1 



-1 1 



1 1 



V 



\ 



/ 1 

1 

V 



These are just the well-known Toeplitz operators D± = T(l =p z) and their transposes are 
denoted by . We also need the infinite diagonal matrix R = diag(|, 1, 1, . . . ) . 



Theorem 5 For sequences of numbers {ctn}.^ 
satisfying 
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define 



A — {aj_k - «j+fc+2)^fc=o - 

Then the following holds true: 
(1) If al = 2afe - ak-i - a^+i , then D+AD^ = . 
(^^y) // = 2afc + afc_i + au+i , then D_AD'^ = A~ . 
(3) If a* - 2a+ + a+_i + a+.i , t/ien D.A+D'^ = 

// a# = 2a^ - a^_i - a^^, , then D+A-Dl = RA*R. 



^j-k + ^j+fc+l)j^fc=0 



# 



(6.5) 



Moreover, if we define a sequence {bn}'^=o and an infinite Hankel matrix by 



1 " 



# 

n-2k 



k=0 



B - ibj+k)Tk=o^ 



(6.6) 



then 



B ^ S^RA*RSl with S: 



# 



/O 



(D 

G) ffl 



(6.7) 



Proof: Before we start with the actual proof, we remark that the various products of the 
infinite matrices make sense in terms of the usual matrix multiplication because the left and 
right factors are always (infinite) band matrices. 

In order to prove the first statement (1) we consider the (j, k) -entries of the following 
(products of) infinite matrices and compute as follows: 



{2aj-k — ttj-k-i — ttj-k+i) — {aj+k+2 — 2aj+fc+i + aj+k) if j,k ^ 1 



[D^ADl] 



[a_k — a_fc+ij — [ak+2 — o-k+i, 
{aj - aj.i) - {aj+2 - aj+i) 

, 0.0 — 0,2 



if J = 0, A; ^ 1 
if j ^ 1, A; = 
if j = A; = 



at 



-,+ 



alk + H+i 

4 + <+i 
On + at 



ifi = 0,A;^l 
ifi ^l,k^O 
iij^k^O 
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Herein we use the fact that 0^ — 0^-1 — 0^+2 + 0^+1 = (20^ — 0^-1- 0^+1) + (2ajt+i — 0^ — 0^+2) j 
a similar identity statement for j , and Oq — ^2 = (2ao ~ 2ai) + (2ai — ag — (22) • Moreover, 
we use the assumption that a„ = a_„ . 

Similarly, we compute the (j, fc) -entry for the product appearing in (2): 



C^Ctj-k + 0,j-k-l + 0,j-k+l) — io,j+k+2 + 2aj+fc+i + a^+fc) 

(a_jt + a_fc+i) — (aj!c4.2 + a^+i) 
(aj + aj_i) - {aj+2 + %+i) 
, ao — 02 



ifj,fc^l 
if j = 0, A; ^ 1 
if j ^ = 
if J = A; = 



- a^+fe+i if j, A; ^ 1 



if J = 0, A; ^ 1 
if J ^ 1,A; = 
if J = /c = 



[^1 



Here we used + afc_i - 0^+2 - a-k+i = C^cik + flfc-i + ak+i) - C^cik+i + Ofc + 0^+2) , a similar 
identity for j , and oq — 02 = (2ao + 2ai) — (2ai + ao + 02) , and again a„ = a_n . 
As for statement (3) we consider 



' (2a+_fc + a+_fe_i + + + 2a++fc + 



[D_A+D^] 



3 

-k+l 



j+k-lJ 



ifj,fc^l 
ifj = 0,fc^l 
if j ^ 1, A; = 
if j = A; = 



U4 



7 On + «( 



ifj,A:^l 
if j = 0, A; ^ 1 
if j ^ 1, A; = 
if j = A; = 



[RA*R] 



Again, we used that = al„ and a* = a*„ . 

Statement (4) can be proven in the same way as statement (3). In fact, if we assume 
all the hypotheses in (l)-(4), then (4) follows with a little algebra from the previous three 
statements (and the from fact that and -D_ commute). 

In order to prove formula (16. 7p first observe that 



# 



00 



(i-j) if i ^ j and i — j even 



otherwise. 



Put tq = 1/2 and r„ = 1 for n ^ 1 . Then the identity B = S^RA^RS'^ can be rephrased 
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as 



m=0 ^ ^ i=0 A;=0 



to hold true for each i, / ^ . These identities are vahd if for each integer s ^ , the 
coefficients for of = ats are the same on the left and right hand side. 

First assume s > . In the right hand side, the coefficient for af = aJtg is equal to the 
sum Ni + N2 + , where 



u I \v 



Ul \V 



fi-2u ^l~2vy 



< M ^ i/2 

O^fe^i i^v ^1/2 

s = j — k s = i — 2u — I + 2v 

N2 = Y ^(hjm.k)r,r,= 

^ M < i/2 

s = k — j s = —i + 2u + I ~ 2v 

N3 = Y ahjm,k)r,rk= Y 

^ i < i ^ u ^ i/2 

s = j + k s = i — 2u + I — 2v 

Therein, we made a change of variables j ^ u = {i — j)/2 and k ^ v = {I — k)/2 . The 
summation is over integer pairs {u, v) . In the above expressions for Ni and N2 we make 
another change of variables v ^ I — v and u ^ i — u to get the expressions 



u I \v 



r i-2u fl-2v 



l\ , V- fi\fl\ 

1^2u-i ^l-2v 



Y^ ( I ( I 'fi~2u f2v-i and 

^-^ \ul \v I ^-^ \UJ \v 

^ M < j/2 \ / \ / i/2 ^ M ^ j 

1/2 1^1 0<i)^Z/2 
s = i — 2u + I — 2v s = i — 2u + I — 2v 

Since there are no indices (m, v) satisfying i/2 ^ u ^ i , 1/2 ^ v ^ I , and s = i — 2u + l — 2v. 
we obtain that Ni + N2 + equals 



(^^(^\=^(i + ls) = l ^^'^ if S < 2 + / and i + / - S even 

^ Vw/ W ~ '^^~\ otherwise. 

^ t; < « 
s = i — 2u + i — 2i) 

This is the desired result since the coefficient for af = dt^ in the left hand side of (16. Sp is 
zero if i + 1 — s is odd and 

2 I I W-s j ~^ \ i+l+s j ) ~ \ i+l-s 
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otherwise. 

In the case s = the coefficient for the term a* on the right hand side of (16.81) equals 
:= N1+N3 = N2+N3 , while it equals ^C,{i+l, 0) on the right hand side. The manipulation 
of the expressions Nk can be done in the same way, with the only difference that in the end 
there are indices which {u, v) satisfying i/2 ^ u ^ i , 1/2 ^ v ^ I , i + I = 2{u + v) . This 
corresponds to a term iV4 , which happen to be equal to . Thus = ^{Ni + ■ • ■ + A^4) 
with Ni + ■ ■ ■ + equaling fl6.9p . This settles the case s = . 

Hence we have shown that the identity fl6.8l) holds, and this implies formula (16. 7p . □ 



In regard to the ffist part of the theorem we remark that the hypotheses in (l)-(4) are 
compatible to each other in the sense that the hypotheses in (1) and (3), as well as those in 
(2) and (4) imply that 

Ofc = 2ak — ak-2 — O'k+2- 

Correspondingly, we have 

DAD^ = RA*R with D = D+D^ = D^D+. 



Elaborating on formula (16. 7p we remark that assuming the hypotheses in (l)-(4), one 

± 

k 



can express the coefficients hn in terms of and as well. We record the corresponding 



results form completeness sake: 



^n = Y.(^^{aUk + 4n-,i-k\ B = S+A'^Sl S^ = S#D. (6.10) 

6„ = ^ (a;_2fe - (hn+i-k)^ B = S-A-S^, = S#D+ (6.11) 

K = ^yA {an~2k - a2n+2-fc), B = SAS^, S = S#D (6.12) 

I.— n ^ j' 



k=0 

The matrices S± and S evaluate as follows: 



/O 

© 
(?) 
Q 
V : 



(D 




Q) 




(I) 







\ 



/ (S) 

-Q 
(?) 
-(?) 
(S 



(I) 
-Q 
(?) 
-(?) 



(?) -ffl 



\ 



■ I 
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s 



/ \ 

(?) - © Q 

(D-(o) (o) 

Q-0 (D-G) Q 
V : -J 

Finally we remark that the recurrence relation (16.61) allows to express the coefficients a* in 
terms of 6„ , 



26n 



a* = a* 



[-] 



Let us now proceed with establishing the identities for the determinants of the finite 
matrices. We restrict to the cases where the symbols are -functions because this is what 
is of interest to us. 

Theorem 6 Let a,a^ G L}{T) he even, and let h E L^[—l,l]. Assume that 



1. a(e 



6(cos^)(2 + 2cose)-i/2(2-2 



cos I 



-1/2 



2. a+(e^^) = 6(cos0)(2 + 2cos0)-V2(2 -2cos0)i/2, 

3. a-(e^^) = 6(cos^)(2 + 2cose)i/2(2-2cose)-i/2^ 

4. a#(e*^) =6(cos0)(2 + 2cos^)i/2(2-2cose)i/2. 

Then, for each 1 , 

det Hn[h] = det (^T„(a) - H^iz-^a)^ = det (^r„(a+) + if„(a 
= det {TJa-) - HJa^^ 



^det (Tn{a*) + H{za*\ 



(6.13) 



Proof. We ffist notice that the hypotheses on the coefficients stated in (l)-(4) of Theorem 
[5] can be rephrased in terms of the corresponding generating functions (see (16.11) and (16.31) ) 
as follows: 



a+(z) = a{z){l- z){l- z-^), 
a*{z) = a+{z){l + z){l + z~^), 



a-{z) = a{z){l + z){l + z-^), 
a#{z) = a-{z){l-z){l-z-'). 



(6.14) 



Here z = e*^ G T . Incidentally, the relations between a , a"*" , a , and a* implied by the 
assumption 1.-4. above are precisely those in (I6.14p . 
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Now assume 4. and compute 



bn = - [ b{x){2xYdx = - f 6(cos0)(2cos0)"sin(^)rf^ 
Ji ^ Jo 

which is precisely the condition (16.61) . 

In order to use the results of Theorem [S] we take the finite sections of the various identities 
(i.e., we consider the n x n upper-left corners of the infinite matrices), 

D+ADl = A+, D_AD^_ = A-, 

D_A+D'^_ = D+A-Dl = RA*R, B = S#RA*RSl, 

and then take the determinants. The crucial point is that , -D_ , and S# are lower- 
trinagular and have ones on their diagonals. The diagonal matrices R give a factor j 
in the determinants. Now it just remains to check that the finite sections of the infinite 
matrices (16. 5p and (16. 6p are indeed the matrices occuring in (I6.13p . But this follows from 
the definitions flO) and (^^. □ 

It is apparent from the proof that if we are only interested in an identity between two types 
of determinants featuring (I6.13p . then it is enough to assume that only the corresponding 
symbols are -functions and that the appropriate relationships between these symbols hold 
(see also (16.141) ). For instance, if we assume a+,a~ G L^(T) and 

a+(^)(l + z){l + z~^) = a~{z){l - z){l - z"^) 

then we can conclude that 

det (T„(a+) + i/„(a+)) = det (T„(a-) - i/„(a-)). (6.15) 

By the way, this relationship between these two type of determinants is not the trivial one 
featuring the "equivalence" between the Cases 2 and 3, which has been pointed out earlier. 

7 Results from the Toeplitz theory 

The idea for this section is that if the a and /5 are any combination of ±|, then we may 
choose an operator of the form T„(a) + Hnib) from our list of identities 1.-4. that has, 
in a certain sense, a nice symbol and find an explicit formula for the determinants of the 
associated Hankel matrices, Hn\b] . This is because for these values of the parameters and 
the right choice of operator, we lose the square root singularities. Fortunately in [B] exact 
formulas for the types of Toeplitz plus Hankel determinants that appear in the previous 
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theorem were found. If we specialize the results to the cases at hand we can state the exact 
formula of the determinants of the matrices Hn\b] . The four different determinants all have 
the form: 

G'[a]"F[a]det(/ + Q„ifQ„), (7.1) 

where F[a\ and G[a\ are certain constants that depend on our choice of parameters for a 
and (3. The last operator determinant involves orthogonal projections Qn = I — Pn , where 
the projections P„ acting on , Z+ = {0, 1, . . . , } , are defined by 

P„(ao, ai, . . . ) = (oq, ai, . . . , a„_i, 0, 0, . . . ). 

The operator K , acting on £^(Z) , is a certain (trace class) semi- infinite Hankel operator. 

The precise reference for the result f l7.ip is Proposition 4.1 and the remark afterwards in 
[6]. Propositions 3.1 and 3.3 in [6j also have to be consulted. For sake of clarification we 
remark that our cases 1.-4. correspond to the cases I-IV in [6] as follows: l.=III, 2.=I, 3. =11, 
4.=IV, where in Case 4, the operators differ by a constant. 

In our case the symbol is (up to a constant) a(e*^) = e~*'^°^^ whence ip = a!!J!^(e*^)a+(e*^) = 
giisme ^ -y^^j^icii occurs in the definition of the operator K . The Fourier coefficients ipk ( ^ ^ ) 
are precisely equal to the value of the Bessel function Jfc(t) of order j with the argument 
t . The precise description of K is as follows: 

1. Let a = = |, then K = Ki where Ki has (j, /c) -entry —Jj^k+2it) . 

2. Let a = —\, 13 = |, then K = K2 where K2 has (j,A;) -entry Jjj^^k+iit) . 

3. Let a = = — |, then K = K3 where has (j,A;) -entry —Jj^k+iit) . 

4. Let a = —\, [3 = — |, then = where has (j,A;) -entry Jj^kif) . 

Here j, k >0. It is known that the operator K is trace class. This is not hard to see since 
for fixed t the entries in the Hankel matrix tend to zero very rapidly. We state the four 
cases below. In all cases our function a in the previous identities is e~^^°^^ times a factor 
of a power of 2 . 

Theorem 7 Let b{x) = (1 - a;)"(l + a;)^e-*^ . 

1. Let a = j3 = J, then 

det Hn[b] = 2-"e*'/^ det (/ + QnKiQn) • 

2. Let a = = |, then 

det H^b] = e*'/*-*/2 det (/ + Qni^aQn) • 
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3. Let a — J, P = — |, then 

det Hn[b] = e*'/«+*/2 det (7 + QnK^Qn) . 

4. Let a — — |, (3 — — |, then 

det i7„[6] = 2"-ie*'/8 det (/ + QnK^Qn) ■ 

This does not quite give us the identity of the original -D„(t) since the above Hankel was 
defined with some extra constants of tt and 2 . So first we adjust for these to yield the 
following. 

Theorem 8 Let b{x) = (1 - + a;)'^e-*^ . 

1. Let a = = |, then 

Dn{t) 

2. Let a — —■^,(3 — |, then 

Dn{t) 

3. Let a — \,(5 — — |, then 

Dn{t) 

4. Let a = — = — |, then 

Dr,{t) = 2-"("-i)-i(27r)"e*'/« det (7 + g„7r4Qn) ■ 

Since Q„ tends to zero strongly and the operator K is trace class, the term det (7 + QnKQn) 
tends to one and the asymptotics are given by the previous factors in each case of the above 
result. 

More precisely, we obtain the following result. 

Theorem 9 Let h{x) = (1 - xY{l + xj^e"*^ . 

1. Let a — ^, (3 — ^, then 

~2-"("+^)(27r)V/l 

2. Let a — — = |, then 

L'„(i)~2-"'(27r)V'/8-*/2. 



= 2-"("+i)(27r)"e*'/« det (7 + QnKiQn) . 
= 2-"'(27r)"e*'/^-*/2 det (7 + Q^K^Qn) • 
= 2-"'(27r)"e*'/«+*/2 det (7 + QnKsQn) . 
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3. Let a = \, (3 = — |, then 

D^„(t)~2-"'(27r)V'/«+*/l 

^. Let a = — = — |, then 

Dnit) ~ 2~"("-i)^i(27r)V'/l 

If we expand det(/ + QnKiQn) using the fact that log det(/ + A) = ti log(/ + A) using just 
the first couple of terms it seems reasonable to conjecture that, for example, 

Dn{t) ~ 2-"(''+^)(27r)V /^e ^2^+3) 

Similar conjectures can be made in the other cases. 

Before ending this section, we conjecture, with the aid the linear statistics formula in 
[T6] and [3] obtained through the heuristic Coulomb Fluid approach [H], that for "general" 
values a and P and for large n 

where 

Here G{z) is the Barnes G -function [2]. 

Finally, we have as a consequence of the previous sections the following remark: Let 
a = |, /3 = |, and let 

(/. = 4log(/ + Q„iri(t/2)Q„). 



Then the function 



t , t. 



(hit) = h -Pi(n, -) 

^ ^ 16 2^^ 2^ 



and thus also satisfies a related Painleve equation. A similar expression can be obtained for 
the other three cases. This once again demonstrates that the most fundamental quantity in 
the theory is the coefficient Pi{n,t). 
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